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Using an equation of motion for a self-gravitating filament, we show 
how galactic spiral arms might be created and sustained. We find that 
the combination of differential rotation of the galactic disk and the self- 
gravity of the arm (as given by the equation) leads to a rotating spiral 
structure. Moreover, using this analysis, we then find a second differential 

^— h equation that explicitly relates this spiral structure to the rotation curve 

of the galaxy - it connects several factors, including spiral shape and pat- 

^ tern speed. We also describe a simple many-body numerical experiment 

that supports our approach. The findings are with consistent with obser- 
vational evidence concerning arm structure and rotation curves, including 
o | leading arm structures. 

One school of thought holds that spiral arms of galaxies cannot be material - 
rotation curves imply that material arms, that is, arms populated with the same 
stars and particles, would wrap up as the galaxy turns [1,2]. Therefore it seems that 
the spiral arm structure must be some kind of density wave - the stars rotate at the 
speeds given by the rotation curve, but somehow cluster along the spiral arms in 
the process. Various mechanisms have been offered for the dynamics of arms [1,3]. 
Others believe that the arms are more transient in nature - breaking up before the 
differential rotation has had a chance to fully wrap them up [4,5,6,7]. Some think 
that a central bar is an important factor in the formation of the arms [8], others 
have put forward the notion that the arms are the result of interactions with another 
galaxy [9]. In any case, the formation and sustenance of spiral arms might fairly be 
called a topic at least open for discussion. Here we offer a simple mechanism that 
can support either transient or long-term arms. 
^ In [10] is introduced an equation of motion for a self-interacting filament, it applies 

in the gravitational case [1 1] . It is 

a =C[2p'T + p K N], 

where a(s) is the filament parametrized by arclength s, p(s) is the mass density 
along a, p'(s) the rate of change of the density with respect to s, T(s) the unit 
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tangent vector, N(s) the unit normal vector and k(s) the curvature along a(s), C is 
a constant, and a dot signifies the time derivative of the motion. 

The derivation of this equation is akin to the Biot-Savart cut-off approach to 
the equation of motion for a line vortex. The idea is that since a one-dimensional 
distribution of mass (a mass density on a space curve) is a divergent object under 
self-gravity, we can analyze the singularity to arrive at the equation. This is a local 
theory, in that it does not account for mass that is close to the test point but which 
is not on the curve, nor for mass that is distal. 

We apply this equation to a spiral arm, and for the moment focus on the k term - 
for momentary illustrative simplicity assume that p' — along the arm. (One could 
assume the arm mass density is given by the exponential decline with radius usually 
assumed for disk galaxies - we do in models below - and note then that since the 
spiral arm goes around the center as well as away, its exponential decline in mass is 
slower than along a ray of the disk, so p' is very small). 

A spiral has decreasing curvature as we follow the arc away from the center. If we 
assume a mass distribution along a spiral, then by the equation of motion we have, 
due to the self-gravitation of the arm, an acceleration at each point in the normal 
direction proportional to the (decreasing) curvature. For a trailing arm spiral, these 
vectors point "in" (have a component toward the center), and backward - against 
the direction of rotation. 

A typical flat or decreasing rotation curve implies shear - differential rotation. In 
order for the spiral pattern to be sustained, each point in the arm must rotate with 
a speed roughly proportional to its distance from the center, for some constant pro- 
portion over the entire arm. The mechanism we propose is that the self-gravitation 
directs the masses in and backwards, more strongly near the center and diminishing 
as we move away, and this compensates for the shear, allowing for a rigidly rotating 
pattern. So, if we are thinking of the arm as a density wave (as opposed to material) 
the orbits of the stars can be thought of as roughly circular with speed given by the 
rotation curve for motion outside the arm. As they enter the spiral arm, they are 
under the additional influence of the arm, directed backward and toward the center. 
Their rotational speed is diminished, creating the density wave that is the spiral arm. 
They then exit the arm to continue around the center. See figures 1 and 2. 

This mechanism is consistent with several aspects of spiral arms. Since the argu- 
ment is largely local, it can support any number of arms, including so-called flocculent 
structures, where there are many vaguely defined arms. There is evidence of a corre- 
lation between spiral pitch angle and shear rate - greater shear tends to give a lesser 
pitch angle [6,7,8]. This is natural in our mechanism - the lesser pitch angle turns the 
normal vector more towards the center, allowing the resonance to be accomplished 
with less force from the arm. See figure 2. 
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There is also evidence of leading (as opposed to trailing) spiral arms - these are 
arms that open into to direction of rotation [9,10]. We can find a similar resonance 
pattern for this sort of arm, though in this case we require a rotation curve such that 
outside layers move past inside ones in the direction of rotation. The observational 
evidence is that leading arms have such rotation curves. In this case our mechanism 
gives a spiral pattern that rotates faster at each radial distance than the speed given 
by the rotation curve. See figure 2. 

Under this mechanism, the spiral arm is a stable structure - an attractor. If the 
self-gravity of the arm, as given by the curvature, is not great enough, then the arm 
winds up further, thus creating greater curvature. On the other hand, if the curvature 
is too great, it defeats the winding, flattening the arm - decreasing curvature. 

We can make a simplifying assumption that allows us to mathematically model 
spiral arm structure and motion in two ways. We run an elementary numerical 
many-body experiment, and we find a differential equation, which depends on the 
rotation curve of the galaxy and the pattern speed, the solution of which is the 
spiral arm (we solve it numerically). The simplifying assumption is that the masses 
keep (approximately) their radius with respect to the center of the galaxy, and so 
the gravitational forces of the arm serve to "locally" speed up or slow them down 
along these roughly circular orbits. We also offer a caveat. In the equation we are 
combining two different sorts of quantities: the gravitational force of the arm as given 
by the filament equation of motion, and rotational velocities. (Of course combining 
forces is natural, as is combining velocities, but not necessarily a mix of the two). 
We argue that because of the force from the arm, a mass passing through the arm 
spends more time in the vicinity of the arm location than it would otherwise, thus 
changing the rotational velocity, and the magnitude of this change is proportional 
to the force. Some justification for this view can be found in figure 3, where we 
have numerically found the trajectory for a mass approaching and passing through 
a filamentary arm. We can see clearly that the amount of time spent near the arm 
depends on the curvature of the arm. We also note that some models treat the 
masses in the galactic disk as if there is a kind of friction among them, that collisions 
and near passes and gas and particle clouds tend to keep masses at the same radius 
moving at the same speed. This tendency can be thought of as a force, pushing 
velocities toward certain values. If viewed this way, we can think of our argument as 
a sort of balance of forces between this 'force' and the gravitational force of the arm. 

The numerical experiment is very simple. We start a group of masses in some 
sort of roughly filamentary distribution, with one end at the center. We then ap- 
ply differential rotation to this distribution, moving the masses in circular orbits, 
according to the rotation speed associated with their radius. We observe that the 
arm is stretched by this, and eventually the masses are so dispersed that there is no 
evidence of the original filamentary structure. We then run the experiment again, 
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this time adding a discretized version of the curvature force (a local action between 
nearest neighbors), at each time step in the experiment. In this case the result is 
that we see persistent spiral structures. For some parameter values these structures 
appear permanent, rotating for as long as the computer runs, for others they seem be 
in the constant process of making and unmaking, finishing several revolutions before 
breaking up then reforming. See figure 4. 

The equation with solutions giving shape and motion of the spiral is the following: 

V(r) + A\\Proj rot ( P KN)\\ + B\\Proj rot (2p'T)\\ = Q p r 

Here r is the radius from the center of the galaxy. V(r) is the rotation curve for 
the galaxy. Q p is the constant pattern speed - assuming the arm structure rotates 
rigidly about the center of the galaxy. If a changing pattern is desired, we can replace 
VL p r with another function of r, Q p (r) . k is the curvature, p the effective density 
at a given station along the arm, p' the rate of change of the density along the arm 
at a given station, Proj rot means that we project this vector onto the tangent of 
the rotation at the given station (recall that our simplifying assumption is that the 
gravitational force of the arm increases or decreases the rotational speed). A, B 
are constants, which may be required because the velocities and the density have 
differing units, curvature depends on scale, etc. The sign of A depends on whether 
or not the arm curve opens toward the direction of rotation - for example for a 
trailing arm A < if the motion given by V(r) is in the counterclockwise (positive) 
direction, so the force from this term is in the clockwise (negative) direction. The 
sign of B depends on whether the density increases or decreases along the curve in 
the direction of the rotation. For a typical trailing arm spiral we have B > 0. 

If we can write all quantities as functions of radius r, then our equation becomes 
an ordinary differential equation with independent variable 9 and dependent variable 
r, so with polar solution r{9). 

Let f3 be the angle of inclination of the arm at a given station (a circle about the 
center of the galaxy has (3 — 0, a ray from the center of the galaxy has j3 — f). 
Then sin((3), cos((3) will project the force, as given by pnN and 2p'T, onto the 
tangent of the rotation. But f3(r) = arctan(^j), and «(r) = - ^ r ~ r 3 r • ■ One way 

of approaching the density is to let the density of the arm simply be the density 
of the disk, that is, we assume that the density of the disk is given by a radially 
symmetric function, so p = Pdisk{r)] this function is usually assumed to be some sort 
of exponential: Pdiskij') = C\e C2T for the portion of the disk not in the galactic bulge 
near the center. Whatever the function is, we would have p\r) = p' disk (r)sin({3), 
because the rate of change of the density along the curve would depend on the angle 
of inclination. 

Then the equation becomes: 
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We also observe that an arm with constant density has p' = (this might be a 
reasonable approximation if the exponential rate c 2 is small enough), so drops that 
term and we have: 
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V(r) has been estimated for many galaxies, as has p(r), at least implicitly, if we 
assume the density in the arm varies radially as the density of the galactic disk does. 
A and Q p are then parameters that, together with initial conditions r(0) and r'(0), 
can be used to fit the arm. That is, we can numerically solve this equation, adjusting 
these values until we have a reasonable fit to images. We have done so for a couple 
well known examples, Messier 101 (a grand design trailing arm spiral - figure 5) and 
NGC 5422 (which is believed to have leading outer arms - figure 6). For Messier 101 
we assumed a flat rotation curve (arbitrarily chosen as V(r) = 2). NGC 5422 has an 
increasing rotation curve - we fit the available data points with a cubic, then used 
that for the V(r). 

Numerical solutions to this equation also reveal the connections between rotation 
curves, pattern speeds, and arm shape in this model. In figure 7 we have used a flat 
rotation curve, fixed A, r(0), and r'(0), and varied the pattern speed Q p . We can 
see that if we knew each of the fixed values, we could, in principle, determine the 
pattern speed from the shape of the arm. 

In this analysis there is nothing surprising about leading arm structures - they are 
simply different values of the parameters and different rotation curves. 

A significant percentage of spiral galaxies are barred spirals - they have a bar in 
the center, and the spiral arms originate at the end of the bar. These bars seem 
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to be reasonably modeled as ellipsoids. We observe that if think of the mass in an 
equidistributed ellipsoid as projected normally onto a line along the major axis, this 
gives a quadratic distribution, the distribution that can be shown, via the above 
equation of motion for a filament, to be the one required for rigid rotation. Moreover 
the rotation curves for matter near the center of the disk often seem to rise roughly 
linearly, as one would expect for rigid rotation. So we could perhaps model the entire 
structure of barred spirals with the filament equation of motion. 

It is well known that there does not appear to be enough matter in a typical disk 
galaxy to hold it together, given the rates at which they appear to spin. This has 
led to the theory of dark matter, additional mass which cannot be 'seen', but which 
provides gravitational stability. The distribution of this mass in a typical galaxy is 
a matter of study. While we do not directly address dark matter here, it may be 
that the model has some implications in the discussion, in that we could ask which 
distributions are consistent with the data (rotation curves and spiral structure) and 
the model. 
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Figure Legends 

Figure 1 

There are three curves, A, B and C. B is the future image of A under a flat rotation 
curve. We can see that the differential rotation has stretched the curve and increased 
the wrapping around the center. The arrows originating on B are normal vectors 
proportional to the curvature, as given by the self -gravitation of the curve. The 
curve C is the image of curve B under this force. C is also the future image of the 
curve A under rigid (no shear) rotation. 

Figure 2 

Schematics of arm 'resonance'. Horizontal lines represent circular motion. The 
more vertical slanted lines represent arms. The first diagram is the trailing arm 
situation. Line A is sent to line B by differential rotation. Line B is sent to line 
C by normal motion (decreasing with distance from the center as curvature does). 
The combination sends line A to line C, rigid translation. The second diagram is 
the leading arm situation. Line A is sent to line B by differential rotation, though 
in this case outside rotation is greater than inside rotation. Line B is sent to line C 
by normal motion (decreasing with distance from the center as curvature does). The 
combination sends line A to line C, rigid translation. In this case the arm pattern 
rotational speed is greater than the disk rotation speeds. The third diagram shows 
that greater shear (consider the distance between lines A and B - it is much greater 
at the bottom than it is at the top) can be compensated for by greater pitch angle 
of the arm. 

Figure 3 

Time spent in/near an arm depends on the curvature of the arm. A. Masses sent 
toward arms of different curvature (k = 0, 1, 2) with same speed and distance from 
arm. The force is calculated as an inverse square integral over the filamentary arm, 
with an e cutoff near the point of intersection of the arm and the trajectory of the 
mass (the trajectory is along the horizontal axis). B. The solutions are approximated 
by simple euler approximation. Vertical axis is position, with the intersection point 
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with the arm, horizontal axis is time. Curve with fastest escape is k = 0, slower is 
k — 1, k — 2 might be said to not have escaped (yet). 

Figure 4 

Snapshots of numerical many-body simulations. The masses are restricted to their 
original radius and so move in circular orbits. In the first shots they are driven only 
by the differential rotation given by a flat rotation curve, and we clearly see the 
decoherence of the original filamentary structure by the winding. In the next shots 
the only change in the system is that we have added a small change in the rotational 
velocities by a discretized version of the curvature between nearest neighbors. The 
persistent spiral structure is evident. 



Figure 5 

Modelling a spiral arm of Messier 101. We have used a flat rotation curve. 
Figure 6 

A. Modelling one of the exterior spiral arms (thought to be leading) of NGC 5422. 
B. For the rotation curve we have used a cubic approximation to data points taken 
from reference 17 (the cubic approximation is pictured). 

Figure 7 

Gallery of spiral arms given by numerical solution to arm equation. Here we have 
fixed flat rotation curve and have varied pattern speed. Of course other combinations 
are possible. 
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